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DOPPLER  STATISTICS  OF  OCEAN  VELOCITY  VARIABILITY 

I.  INTRODUCTION 

In  the  past  decade,  the  oceanographic  community  has  shown  a  great  deal  of 
interest  in  the  use  of  high-frequency-Doppler  sonar  to  measure  the  ocean 
velocity  field.  Pinkel  (1979,  1981)  successfully  used  a  87.5-kHz,  3-kW  system 
mounted  on  the  R/P  FLIP  to  measure  the  velocity  of  internal  waves.  Regier 
(1982)  and  Joyce  et  al.  (1982)  demonstrated  that  a  Doppler  profiler  can  be 
employed  underway  from  a  ship  to  measure  current  shear.  L'Hermitte  (1973)  and 
Farmer  and  Crawford  (1983)  used  coherent  Doppler  sonar  techniques  to  infer 
ocean  currents.  Recent  review  articles  in  ocean  Doppler  sonar  techniques 
include  Pinkel  (1979)  and  Mathews  and  Hicks  (1981). 

Despite  all  of  this  recent  work  in  the  area  of  Doppler  sonar,  there  has 
been  only  limited  research  into  the  relationship  between  Doppler  signal 
statistics  and  ocean  velocity  statistics,  and  into  the  intrinsic  limits  of 
performance  of  oceanographic  Doppler  sonar  systems  imposed  by  the  nature  of 
the  velocity  field  itself,  in  particular,  by  the  spatial  statistics  of  the 
velocity  field.  Oceanographic  Doppler  systems  typically  are  limited  to 
estimating  the  Doppler  shift  associated  with  the  component  of  velocity  along 
the  mean  acoustic  ray  path  of  the  beam  pattern,  averaged  over  the  scattering 
volume.  The  effect  of  a  random  component  of  the  velocity  field  and  the  effect 
of  the  component  of  velocity  transverse  to  the  mean  beam  pattern  axis  is  to 
produce  Doppler  spectral  broadening  that  is  usually  treated  as  noise.  (See, 
for  example,  Pinkel  (1981),  who  employed  the  spectral  estimation  techniques  of 
Runnier  (1968)  to  obtain  internal  wave-induced  Doppler  velocities.) 

The  problem  of  Doppler  broadening  by  random  velocity  fluctuations  and  by 
flow  perpendicular  to  a  finite  sized  beam  pattern  has  drawn  some  interest  both 
in  the  general  ultrasound  literature  (see  Mathews  and  Hicks,  1981)  and  in  the 
laser  flow  literature  (see  Edwards  et  al.,  1971).  Doppler  broadening  by  a 
random  velocity  field  has  been  considered  by  a  number  of  authors,  including 
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Green  (1964),  Rasmussen  and  Head  (1978),  Farmer  (1983),  and  L'Hermitte 
(1983).  Their  results  show  that  this  type  of  Doppler  broadening  is 
proportional  to  the  rms  value  of  the  fluctuating  velocity. 

Brown  and  Clifford  (1973),  Brown  (1974),  and  Sullivan  and  Kemp  (1979), 
among  others,  examined  Doppler  broadening  effects  due  the  advection  of  the 
component  of  the  velocity  field  perpendicular  to  the  beam  pattern  axis. 
Newhouse  et  al.  (1980)  and  Albright  (1976)  examined  the  combination  of  both  of 
the  above  effects. 

The  purpose  of  this  report  is  to  develop  a  model  that  will  be  useful  in 
obtaining  the  relationship  between  the  space/time  statistics  of  the  ocean 
velocity  field  and  those  of  the  returned  acoustic  signal  for  a  pulsed,  high- 
frequency,  narrowbeam,  monostatic  sonar  system  —  the  type  of  system  commonly 
employed  in  oceanographic  measurements.  Special  attention  is  given  to  the 
statistics  of  the  pulse-to-pulse  correlation  function,  in  particular,  its 
variance,  which  is  a  fourth-order  moment.  Also,  the  utility  of  using  the 
statistics  of  the  returned  signal  to  infer  information  about  the  spatial 
structure  of  the  microstructure  velocity  field  is  analyzed. 
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II.  SCATTERING  MODEL 

Consider  the  signal  received  by  a  narrowbeam,  pulsed,  monostatie  sonar 
resulting  from  the  acoustic  scattering  of  discrete  particles  embedded  in  a 
randomly  moving  fluid.  The  following  assumptions  are  made: 

1.  The  frequency  is  sufficiently  high  that,  in  the  farfield,  straightline 
geometric  acoustic  propagation  applies. 

2.  The  scatterers  are  assumed  to  follow  the  fluid  motion  perfectly. 

The  returned  signal,  then,  is  a  linear  sum  of  reradiated  replicas  of  the 
transmitted  signal,  modulated  in  amplitude  by  individual  particles  strength, 
transmission  loss,  and  the  transmitter/receiver  beam  pattern,  and  in  phase  by 
the  location  of  each  scatterer.  Only  the  farfield  case  will  be  examined. 

A  monostatic  sonar  transmits  a  series  of  identical  pulses  $j(t)  of 
duration  Tp  repeatable  at  time  interval  Tr.  Each  pulse  is  expressed  as  the 
Fourier  series  over  the  pulse  duration  Tp.  Thus, 

~~  ao 

4>x(t)  »  E  Yk  exp(iufct)  for  -Tp/2  <  t  <  Tp/2  , 
k— « 

(la) 

$T(t)  -  0  for  Tp/2  <  t  <Tr  -  Tp/2  , 

and 

$T[t  -  (n  -  l)Tr]  -  *r(t)  for  0  s  t  -  (n  -  l)Tr  £  Tr  ,  (lb) 

n  -  l,  2,  3,  - 

where  «k  -  2irk/Tp  -  kA«.  The  pulse  is  taken  as  narrowband  about  some  center 
frequency  ojq,  with 


where  r  «  q,  with  the  bandwidth  B  given  by  B  -  2rAo>.  The  contribution  of 
ambiguous  scattering  volumes  in  the  returned  signal  is  ignored.  (This  can  be 
accomplished  by  either  taking 


Tr  >  2RJnax/c  »  Tp  , 

where  R^ax  is  the  maximum  range  at  which  the  returned  scatterer  signal  is 
above  some  prescribed  detection  threshold,  or  by  the  use  of  a  series  of 
narrowband  pulses  with  a  sufficiently  wide  range  of  center  frequencies.) 

The  returned  signal  $(t)  of  a  narrowbeam,  monostatic  sonar  system 
operating  in  the  geometric  acoustics  regime  is  just  the  sum  of  the 
contributions  of  the  individual  discrete  scatterers  and  can  be  represented  by 

4>(t)  -  S  dr  b2(n)a(r)p(r,t)<frr(t  -  2r/c)  , 

-  /  dr  /  b2(a)a(r)p(r, t)  I  ^  exp  i [ufcf t-2r/c) J  ,  (2) 

where  b(n),  n  -  r/r  is  the  beam  pattern  function.  The  summation  over  K  will 
be  understood  to  mean,  subsequently,  over  all  k,  from  +®  to  -®.  The  radial 
integral  in  equation  (2)  is  also  understood  to  be  constrained  by  the  pulse 
duration,  i.e., 

c(t  -  Tp)/2  <  r  <  c(t  +  Tp)/2  . 

The  discrete  scatterer  density  function  p(r,t)  is  defined  as 
N 

p(r,t)  -  E  6[r  -  r(Xj ,  t  -  t j ) ]  ,  (3) 

j-1 


where  r(X,t)  is  the  vector  position  (at  time  t)  of  a  fluid  element.  This 
representation  is  termed  Lagrangian  in  fluid  dynamics  (see  Lamb,  1932).  The 
position  at  initial  time  t  -  0,  namely  X,  is  defined  by  X  -  r(X,0).  For 
scatterers  being  advected  with  the  fluid,  Xj  represents  the  initial  position 
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of  the  jth  scatterer.  The  factor  t  -  tj  in  r  of  equation  (3)  is  the  time  at 
which  acoustic  energy  was  reflected  from  the  jth  scatterer,  which  was  located 
at  the  time  of  scattering  (i.e.,  t  -  tj)  at  location  r(Xj ,  t  -  tj). 

The  discussion  is  confined  to  the  case  where  the  scatterers  follow  the  fluid 
perfectly.  The  quantity  <x(r)  in  equation  (2)  is  taken  as  real  and  represents 
the  signal  amplitude  resulting  from  the  combined  effects  of  transmission  loss 
and  scattering  strength. 

Consider  the  received  signal  ^(t)  due  to  the  mth  pulse.  From  equations 
(2)  and  (lb),  it  follows  that 

♦m(t)  =*  /  dr  E  tk  b2(n)a(r)p( r,  t )  exp{iuk[t-(m-l)Tr-2r/c]}  .  (4) 

V 

The  subscript  V  in  equation  (4)  indicates  that  the  volume  integral  is  over 
scatterers  ensonified  by  the  mth  pulse.  A  range  gate  is  now  applied  to  ^(t) 
and  the  returned  gated  function  ^(t')  is  defined  as 

*m(t' )  -  h( t  )^m( t )  ,  (5) 


where 


t  -  t  -  (m  -  1  )Tr  -  Xu]  , 


and 


h  -  h( t ' )  for  -T[/2  ^t’i  T[/2  , 
h  -  0  otherwise. 

Now,  take 

Ti  «  Tr  . 
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Here,  T,  is  the  duration  of  the  range  gate  and  Tm  is  the  time  delay  of  the 
range  gate.  Substitution  of  equation  (4)  into  (5)  yields 

*m(t)  -  h(t')  J  dr  Z  b2(n)ct(r)p[r,  t '  +  (m  -  l)Tr  + 

V 

exp[i«k(t'+tm-2r/c)]  , 


s  h( t  )<x( Rm)  /  dr  Z  b2(n)p[r,t'  +  (m  -  l)Tr  +  xm]<|rk 
V 

exp{i«k[ t ' +2(Rffl-r ) ]/c}  ,  (6) 

where  Rfl,  -  CTm/2  »  cTp/2  ,  The  signal  amplitude  factor  in  equation  (6)  has 
been  taken  to  be 

<x(r)  3  a ( Rnj )  . 

The  volume  integral  in  (6)  at  any  instant  t'  (see  figure  1)  is  taken  over  the 
range  interval 

(7) 

rt  -  c( t ’  -  Tp/2 )/2  ,  (8a) 


rx  §  r  -  5  r2  , 

wi  th 


and 


r2  -  c(  t  ’  Tp/2  )/2  .  (8b) 

Equation  (6)  represents  the  starting  point  for  Doppler  processing.  Pinkel 
(1981)  estimates  Doppler  shift  by  calculating  an  incoherent  average  of  the 
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Figure  1.  Scattering  Volume  Contribution 
(The  radial  contribution  of  scatterers  for  the  mth  pulse  transmitted  during  the 
time  interval  (m-l)Tr-Tjj/2  £  t  £  (m-l)Tr+Tg/2 ,  where  Tr  is  the  repetition 
rate,  Tp  is  the  pulse  duration,  and  t'  refers  to  the  receiver  gate  time 
defined  over  the  gate  interval  0  £  t’sS  Ti  -  Tp.) 


mean  of  the  spectral  estimate  of  ^(t')  over  the  gated  time  interval  Tj. 

(This  is  a  technique  first  suggested  by  Rummler,  1968.)  The  Doppler  velocity 
sensitivity  of  each  estimate  in  the  average  is  defined  by  the  effective  pulse 
duration  and  is  improved  by  performing  an  incoherent  average  over  spectral 
estimates  of  a  set  of  pulses  over  a  time  period  such  that  the  estimated 
Doppler-shifted  velocity  remains  approximately  constant. 
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In  this  report,  the  pulse-to-pulse  statistics  of  4*m( t ' )  will  be 
examined.  The  statistics  of  the  problem  enter  through  the  Lagrangian  position 
vector  r(Xj.t)  in  the  scatterer  density  function  p(r,t)  given  by  equation 
(3).  One  can  decompose  r(Xj,t)  as 

r(Xj,t)  -  Xj  +  dj  ,  (9) 

where  dj  =*  d(Xj,t)  is  the  relative  displacement  of  a  scatterer  from  its 
initial  position  Xj .  Note  the  difference  in  the  statistics  of  Xj  and  dj . 

Xj  refers  to  the  position  of  scatterers  at  a  particular  instant  of  time  (in 
this  case,  the  initial  instant  of  time  t  -  0),  while  dj  is  the  displacement 
from  that  position  due  to  advection  by  the  fluid  velocity  field.  It  will  be 
assumed  that  Xj  has  an  equiprobability  of  initially  being  anywhere  in  the 
spatial  domain  of  the  problem  (or,  equivalently,  each  of  the  N  scatterers  are 
placed  equiprobably  throughout  the  volume).  Note  that  the  statistics  of  Xj 
and  dj  are  independent  and  that,  therefore,  for  any  function  F(Xj,dj),  the 
ensemble  average  <F(Xj,dj)>  can  be  written  as 

<F(Xj,dj)>  -  <<F>x>d  -  <<F>d>X  - 

where 

<F>X  -  X  dXj  FP(Xj)  , 
and 

<F>d  -  /  ddj  FP(dj)  , 

respectively,  and  P(  )  indicates  the  probability  density  function.  The  time 
average  is  also  defined  as 

T 

{(•••))T  -  (1/T)  X  dt  (•••)  , 

0 
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and  the  space  average  as 


L 

[(  -)]L  -  (1/L)  /  dx  (•••)  . 

0 

It  is  straightforward  to  show  from  equation  (6)  that  the  time  average  of  the 
received  signal,  the  first-order  time  statistic,  is  zero: 

[*kn(t  ’  JIij  -  0  .  (10) 

For  the  second-order  statistics,  the  pulse-to-pulse  gated  temporal  correlation 
function  S  is  defined  in  the  usual  fashion  and  can  be  written  as 

S'nm  -  [*n(t,)«to(t,)]Tl  *  a(Rn)a(Rm)<h2(r  )  f  dr’  f  dr"b2(n’  )b2(n" ) 

Vn  Vm 

E  W'l  expli^k-MOt  ’+aifuiCRn-r’ )-aii(Rn,-r")]/c} 
k,  1 

p[r  >t  +(n-l)Tr+Tn]  p[r  ,t  +(m-l)Tr+Tm]>Tj  ,  (11) 

where  the  scattering  volumes  Vn,  Vn  from  equation  (7)  limit  the  range  of 
integration  of  (11)  to 

rl  -  r  -  Rn  *  ?2  * 


and 


ri  S  r”  -  Rn  S  r2  , 

respectively.  Note  from  equations  (8)  that  r^  -  c(t'  -  Tp/2)/2  and  r2  - 
c ( t '  +  Tp/2 )/2 .  In  the  appendix,  a  detailed  derivation  is  carried  out  for  the 
normalized  pulse-to-pulse  correlation  function  S^,  -  S'nn/S'nn  for  the  case  in 
which  the  largest  scale  of  the  fluctuation  field  is  smaller  than  the  range 
dimension  of  the  scattering  volume.  Generalization  of  these  results  to  the 
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case  where  the  scattering  volume  is  smaller  than  the  largest  scale  of  the 
microstructure  field,  as  discussed  in  the  appendix,  is  straightforward  but 
algebraically  complicated.  The  final  result  of  the  appendix  derivation  is 
equation  (A-31);  namely, 

Snm  -  z  tit* l  S  dnjjb4 ( njj )  exp[ i2wivh(m-n)Tr-nh/c] 

1 

exp{io)i[2vz(m-n)Tr/c-AT]}  S  d5g(£)  exp{iw[2v ' z(m-n)Tr/c]}  ,  (12) 

where 

v'z  *  v’z(z)  - 
and 

5  *  (z  -  R)/Az  , 

with  Az  -  cTp/2  and  R  -  Rn  .  Also,  vz  is  the  mean  velocity  averaged  over  the 
scattering  volume,  equation  (A-26),  in  the  direction  of  the  mean  ray  axis  of 
the  beam  pattern;  v'z  is  the  random  velocity  in  the  direction  of  the  mean  ray 
axis  of  the  transducer,  and  vjj  is  the  vector  component  of  velocity 
perpendicular  to  the  mean  ray  axis  of  the  beam  pattern.  The  range  gate  used 
in  deriving  equation  (12)  takes  Tj  -  Tp 

and 

h  -  1  for  -T{/2  St'S  T\/2  , 
h  -  0  otherwise. 

It  should  be  noted  that  equation  (12)  has  appeared  in  a  variety  of  forms  in 
the  scientific  literature.  See,  for  example,  Green  (1964),  Edwards  et  al. 
(1971),  Albright  (1976),  and  Farmer  (1983). 
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Equation  (12)  is  further  simplified  by  writing  it  in  the  form  of  the 
product  of  three  terms,  namely, 


Snm  ”  E  =iK iQl  , 


(13) 


where 


=1  -  tl+'l  (14) 

(with  normalized,  namely,  Z  =i  -  1)  is  the  1th  spectral  component  of  the 
transmitted  pulse;  and  where 

Ki  -  /  dnhb4(nh)  exp[ i2«ivh(m-n)Tr*nh/c]  (15) 

represents  the  contribution  from  the  component  of  velocity  perpendicular  to 
the  mean  ray  axis;  and 

Ql  -  exp(i2cji(m-n)Tr(v2-cAt/2)/c]  /  d£g(£)  exp[ i2<aj (m-n)Trv’z/c]  (16) 

is  the  contribution  to  the  correlation  function  from  the  component  of  velocity 
along  the  mean  ray  axis  of  the  transducer,  including  both  mean  and  random 
components.  Note  that,  if  v'z  -  0  and  v^  «  0,  substitution  of  equation  (16) 
into  (13)  yields 

Snm  -  1  'm'lQl  -  E  *1**1  exp{i  exp[i2ui(m-n)Tr(vz-cA‘c/2)/c]}  •  U7) 

Equation  (17)  is  the  theoretical  basis  for  a  coherent  Doppler  estimate  of  vz. 


11/12 
Reverse  Blank 


III.  INPUT  OF  FLUID  VELOCITY  STATISTICS 


Consider  the  role  of  the  mean  along-beam  velocity  component.  From  the 
form  of  equations  (13)  and  (16),  Sj^  has  a  maximum  if 

A-c  -  2vz(m  -  n)Tr/c  ,  (18) 

which  can  be  interpreted  as  an  along-beam  alignment  of  the  mth  and  nth 
scattering  volumes,  which  are  separated  by  the  time  interval  (m  -  n)Tr. 
Equivalent ly,  the  phase  difference  of  the  two  received  signals  produced  by  vz 
in  equation  (16)  can  be  compensated  for  by  adjusting  At  to  maximize  Sun,.  The 
use  of  equation  (18)  represents  a  coherent  Doppler  estimate  of  ihe  Lean 
velocity  in  the  direction  of  the  sonar  axis  (L'Hermitte,  1983).  The  use  of 
equation  (18)  in  (16)  results  in 

Ql  -  S  d^gU)  exp[i2ui(m-n)Trv’z/c]  .  (19) 

Consider  the  factor  Kj  in  equation  (13),  which  from  equation  (A-10)  is 
recognized  as  a  fourfold  convolution  integral  of  the  aperture  function  w: 

K1  ”  w[vh(m  -  n)Tr  -  x']*w*w*w  ,  (20) 

where  the  convolution  operator  is  defined  as 

U(x  -  x’ )*V  -  /  dx'U(x  -  x’)V(x’). 

It  is  straightforward  to  show  that  if  the  weight  w  has  characteristic  lengths 
Lx,  Ly,  then  the  fourfold  convolution  of  w,  i.e.,  equation  (20),  also  has 
characteristic  lengths  L*,  Ly.  Thus  Kj,  which  is  the  sole  contributor  to  the 
pulse-to-pulse  correlation  function  in  the  absence  of  any  radial  velocity 
fluctuations  on  the  scale  of  the  scattering  volume  and  smaller,  has  a 
characteristic  length  scale  of  the  dimensions  of  the  sonar  aperture. 

Therefore,  K[  (and,  hence,  Snuj)  becomes  small  when 

vh(m  -  n)Tr  >  L  ,  (21) 
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where  L  -  (Lx,  Ly)  is  the  characteristic  vector  length  scale  of  the  sonar 
aperture.  Thus,  v^,L  can  be  combined  to  form  a  characteristic  time  of 
decorrelation. 

Note  from  equation  (Id)  that  Qi  is  a  spatial  average.  Consider  the 
circumstances  under  which  this  spatial  average  can  be  used  in  place  of  an 
ensemble  average  over  displacement  statistics,  i.e,  when 

Ql  -  <Ql>d  •  (22) 

Assuming  Gaussian  statistics  and  a  homogeneous  turbulence  (Hinze,  1959)  model 
for  the  microstructure  velocity  field  yields  from  equation  (19) 

< Q 1  >d  ”  X  dCg(5)<exp[i2ui(m-n)Trv,z/c]>tj  -  exp  -x2  ,  (23) 


where 


x  **  (m  -  n)Tr/xo  . 


(24a) 


with 


(to)"2  *  2iq2  <(v'z)2>d  ,  (24b) 

where  iq  -  a>i/c.  Note  that  xq  represents  the  characteristic  decorrelation 
time  associated  with  vz  —  the  radial  velocity  fluctuations  on  the  scale  of 
the  scattering  volume  and  smaller.  It  should  be  noted  that  results  similar  to 
equation  (23)  have  been  derived  by  a  number  of  authors;  see,  among  others, 

Green  (1964),  Edwards  et  al.  (1971),  Albright  (1976),  and  Farmer  (1983). 

To  examine  the  validity  of  equating  the  spatial  average  with  the  ensemble 
average,  i.e.,  equation  (22),  consider  the  variance  of  Qi  defined  by 


9  -  < (Qi  -  <Ql>d)2>d 

-  <Q l >d2  SS  dCdC’g(C)g(C' ){exp[2x2Rv(z-Z' ))-l)  , 


(25) 


where  £  is  defined  in  equation  (12)  and  Rv  is  the  normalized  fluctuating 
z-component  velocity  correlation  function: 


Rv(z  -  z)  -  <vz(z)v'z(z' )>d/<(v'z)2>d  .  (26) 

The  spatial  dependence  of  Rv  on  the  difference  variable  z  -  z’  results  from 
the  assumption  of  homogeneity  for  the  fluctuating  field  vz.  Note  that  9, 
equation  (25),  is  a  fourth-order  statistic  of  the  received  signal. 

By  definition,  |RV|  S  1,  with  the  limiting  value  Rv  -  0  corresponding  to 
uncorrelated  variability  on  the  scale  of  the  scattering  volume  and  smaller  for 
which  case,  from  equation  (25),  9  -  0  and 


Q  -  <Q>d  • 

However,  in  general,  Rv  *  0  and,  thus,  0*0,  and  Q  *  <Q>d  • 


As  a  simple  example  illustrating  the  nature  of  equation  (25),  set  g(£)  - 
1  (even  though  this  is  not  true  and  g  is  given  by  the  triangle  function,  the 
difference  in  results  being  a  constant  on  the  order  of  l)  and  let  Rv  be  given 
by  the  simple  expression 


Rv(z)  -  1  -  z/1  for  z  £  1  , 
Rv(z)  -  0  for  z  >  1  , 


(27) 


where  1  is  interpreted  as  the  characteristic  scale  of  the  variability.  This 
simplification  still  retains  the  essential  physics  of  the  velocity  correlation 
function  and  its  role  in  equation  (25).  The  resulting  9  is  plotted  in  figure 
2  versus  x  for  different  values  of  the  nondimensional  parameter  xq  -  (1/Az), 
the  ratio  of  the  velocity  variability  length  scale  to  the  along-range 
dimension  of  the  scattering  volume.  Figure  2  also  shows  a  plot  of  (<Ql>dJ2  - 
fluid  velocity  field  in  the  z-direction  becomes  very  small  compared  with  the 
exp  -2x2  versus  x.  Note  that  as  xq  **  0,  i.e.,  the  spatial  length  scale  of  the 

15 


Figure  2.  Second-  and  Fourth-Order  Moments 
(The  solid  line  is  the  function  <Ql >d *  equation  (23);  the  symbols  refer  to 
families  of  0,  equation  (25),  for  different  ratios  of  xq  -  1/Az,  with  g(£)  - 
1  and  equation  (27)  as  the  model  for  microstructure  correlation  function.) 


scattering  volume  dimension  Az,  0-0.  This  implies  that  the  variability 
characterized  by  scales  that  are  small  relative  to  the  scattering  volume  (in 
the  z-direction)  result  in  small  0  and,  hence,  Q  s  <Q>jj.  As  xq  "  l,  0 
increases  and  Q  *  <Q>d*  From  the  general  form  of  equation  (25),  note  that 
different  spatial  scales  of  the  velocity  field  contribute  to  0  for  different 
x,  the  largest  scales  contributing  relatively  more  for  small  x  and  the 
smallest  scales  contributing  more  for  larger  x.  To  understand  this,  note  that 
from  equation  (25) 


I 
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f dO/d(x2 ) Jx-0  “  2X  , 


(28) 


where  X  is  the  integral  scale  defined  by 
x  -  SS  dCdi-'gU)gU’)RvU  -  z  )  . 

The  integral  scale  is  a  measure  of  the  largest  scale  of  the  variability.  (See 
Hinze  (1959)  for  a  discussion  of  integral  scale  in  turbulence  theory.)  This 
behavior  is  indicated  through  the  family  of  different  xq  values  in  figure  2  by 
noting  the  large  change  in  the  slope  of  0  near  x  -  0,  for  increasing  xq. 

At  large  times  x,  only  the  smallest  values  of  Rv(z  -  z)  (and,  hence, 
smallest  scales)  contribute,  since  Rv  is,  in  general,  a  monotonical ly 
decreasing  function  for  realistic  ocean  microstructure  fields.  Using  the 
lowest  order  expansion  of  Rv(z  -  z),  namely, 

Rv(z  -  z' )  -  1  -  |z  -  z-|2/(21v2), 

where  lv  is  termed  the  turbulent  microscale,  which  is  a  measure  of  the 
smallest  scale  of  the  variability  field  (see  Hinze,  1959,  or  Tatarski ,  1971), 
in  equation  (25)  yields  for  x  >  1  the  approximation 

<(Ql)2>d  -  SS  d5drgtt)g(C)  exp[-(z-z’)2^2/lv2]  *  iv/(^z)  .  (2g) 

Note  that  while  the  function  Qj  contains  information  on  the  variance  of  the 
velocity  field,  9  or  <(Q])2>d  contains  additional  information  on  the  spatial 
structure  of  the  velocity  field.  Conversely,  estimates  of  both  Qi  and  0  from 
measurements  might  allow  estimates  not  only  of  the  variance  <(v'z)2>(j  but  also 
of  the  velocity  correlation  function  Rv  and,  in  turn,  its  Fourier  transform, 
the  velocity  wavenumber  spectrum. 
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IV.  OCEAN  APPLICATION 


The  theory  developed  in  sections  II  and  III  is  applied  here  to  the  case 
of  ocean  velocity  statistics  modeled  as  microstructure.  The  notation  and 
coordinate  system  of  the  earlier  sections  are  continued.  Note  that  the 
z-coordinate  is  along  the  mean  ray  axis  of  the  beam  pattern.  The  dominant 
velocity  field  transverse  to  the  beam  axis  is  in  the  horizontal.  As  a 
concrete  example,  consider  a  moored,  pulsed,  monostatic,  sonar  system.  A 
comparison  will  be  made  between  the  operational  parameters  for  a  Doppler 
velocity  estimate  and  a  velocity  estimate  that  would  be  inferred  through  the 
use  of  the  higher  order  statistic  < (Q 1  )2>d »  give*1  through  equation  (25),  along 
with  the  second-order  statistic  <Ql >d*  ^rom  equation  (23).  This  latter 
estimation  technique  is  called  stochastic  Doppler  because  it  involves  the 
direct  relationship  between  the  ensemble  statistics  of  the  acoustic  field  and 
the  ensemble  statistics  of  the  oceanographic  velocity  field,  in  contrast  to 
traditional  Doppler,  which  estimates  velocity  in  real  time. 

The  inverse  relationship  between  the  statistics  of  the  velocity  field  and 
0  or  <(Qi)2>d  is  not  developed  in  this  section,  but  is  left  for  a  subsequent 
report.  Rather,  the  limits  of  performance  of  the  both  traditional  Doppler  and 
this  new  stochastic  Doppler  are  examined  in  terms  of  the  minimum  spatial 
resolution,  maximum  velocity  sensitivity,  and  ranges  potentially  achievable, 
given  system  noise  and  the  variance  of  an  estimate  caused  by  the  finite  number 
of  degrees  of  freedom  intrinsic  to  the  measurement  of  microstructure. 

In  examining  the  constraints  imposed  on  real-world  operation,  other 
factors  will  be  ignored  for  simplicity.  These  include  the  individual 
scatterer  strength  and  intrinsic  phase  variability  (apart  from  relative 
location,  which  is  considered  in  this  work),  scatterer  motion  relative  to  the 
fluid,  and  phase  fluctuations  from  propagation  through  random  index-of- 
refraction  variability.  (Note  that  it  may  be  possible  to  infer  scatterer 
motion  relative  to  the  fluid  to  the  resolution/sensitivity  discussed  below  if 
such  motion  is  statistically  different  from  the  fluid  motion.)  In  addition, 
the  reader  is  reminded  that  an  incoherent  scatterer  model  has  been  used  with  a 
large  number  of  scatterers  per  unit  scattering  volume. 
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Small-scale  ocean  variability,  which  includes  internal  waves, 
f inestructure ,  and  microstructure,  is  here  defined  as  variability  on  the  order 
of  and  less  than  internal  waves,  i.e.,  spatial  scales  of  tens  of  kilometers 
horizontally,  hundreds  of  meters  vertically,  time  scales  between  the  inertial 
period  (on  the  order  of  24  hours  at  mid-latitude)  and  the  buoyancy  period 
(typically  10  to  20  minutes).  A  "canonical"  vertical  wavenumber  spectrum  of 
velocity  variability,  based  on  observed  data,  is  shown  in  figure  3  from 


Gargett  et  al.  (1981).  Note  the  regimes  of  internal  waves,  f inestructure,  and 
microstructure.  Implicit  in  this  model  is  that  microstructure  is  isotropic 
with  a  characteristic  (largest)  scale  size  1,  and  a  dissipative  (or  smallest) 
scale  size  lv.  Note  that  1  s  meters  and  lv  s  centimeters  for  ocean 
microstructure  and  that  the  spectral  level  of  velocity  microstructure  varies, 
following  the  three  curves  (I,  II,  III)  of  figure  3. 


Figure  3.  Observed  Vertical  Wavenumber  Spectrum  of  Small-Scale  Ocean  Velocity 
(From  Gargett  et  al.  (1981).  Shown  are  three  levels  of  microstructure, 
labeled  I,  II,  and  III.  Also  indicated  are  regimes  of  internal  waves, 
f inestructure,  and  microstructure.) 
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The  Doppler  approach  is  considered  first.  This  analysis  is  similar  to 
that  of  Pinkel  (1981)  for  the  internal  wave  range  used  in  the  R/P  FLIP'S 
87.5-kHz  sonar  system,  except  that  the  coherent  Doppler  technique  is  used 
(rather  than  the  Rummler  (1968)  technique),  since  the  problem  has  been 
formulated  in  that  manner.  It  can  be  shown  that  the  two  techniques  yield  the 
same  level  of  sensitivity/precision  for  the  same  number  of  independent, 
real-time,  velocity  estimates. 

The  estimate  of  Doppler  velocity  is  obtained  from  equation  (18)  via  the 
shifted  time  delay  At  -  2vz(m  -  n)Tr/c,  where  vz  is  the  spatially  averaged 
mean,  defined  by  equation  (A-26)  in  the  appendix.  Using  equation  (17)  for  the 
case  where  v'z  is  negligible,  which  would  occur  if  the  scattering  volume  were 
sufficiently  small,  then  from  equation  (16)  Qi  -  1  when  At  -  2vz(m  -  n)Tr/c. 
However,  the  contribution  of  Kj  in  equation  (13),  which  contains  the  effect  of 
motion  perpendicular  to  the  mean  beam  axis  and  represents  a  Doppler  spread 
effect,  must  be  considered.  (For  examination  of  this  effect,  see  (among 
others)  Green,  1964,  and  Brown,  1973.)  As  discussed  in  section  III,  Ki  is 
nonzero  over  a  time  scale  L/v*,  where  v*  is  the  characteristic  velocity  of 
flow  perpendicular  to  the  mean  beam  axis,  and  L  is  the  characteristic 
dimension  of  the  sonar  in  the  direction  of  v*.  Note  also  that  if  a  mean 
shear,  q  -  |^hvzl -  °f  the  along-beam  velocity  is  present,  it  is 
straightforward  to  show  that  there  is  a  decorrelation  effect  of  the  same  form 
as  equation  (20)  for  the  transverse  component  of  flow.  The  same  scaling  for  a 
characteristic  decorrelation  time  occurs;  namely,  L/v*,  with  v*  s  qR,  R  being 
the  range.  This  can  also  be  derived  from  equation  (A-25)  using  (A-10)  by 
including  the  shear  term  in  the  expansion  for  the  mean  velocity  in  equation 
(A-26),  whence,  with  this  substitution,  the  shear  term  q  has  the  same  effect 
as  the  mean  component  of  flow  perpendicular  to  the  mean  ray  axis. 

Thus,  for  an  estimate  of  velocity  from  the  pulse-to-pulse  correlation 
function  S^,  the  time  delay  between  pulses  must  be  restricted  to 

t  -  (m  -  n)Tr  <  Tujjm  -  L/v*  ,  (30) 

where  v*  is  taken  as  the  maximum  of  qR  and  the  transverse  velocity  component 
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(perpendicular  to  the  mean  ray  axis),  and  where  equation  (17)  can  be  used  to 
estimate  the  Doppler  time  delay  At.  Note  from  equation  (30)  that  the  maximum 
number  of  correlated  sequential  pulses  p  is  given  by 


p  -  (m  -  n)max  -  L/(Trv* )  . 

The  maximum  time  delay  At  -  ATmax  over  which  an  estimate  of  the  mean  velocity 
vz  can  be  made  is  on  the  order  of 

AT^ax  -  2vzpTr/c  a  2L(vz/v')/c  .  (31) 

To  resolve  this  time  shift,  which  is  estimated  by  maximizing  Snm  (and  hence 
Qj)  given  by  equation  (16),  namely, 

Qj  ■  exp[ i2«i(m-n)Tr(vz-cAT/2)/c]  , 

it  is  necessary  that  the  fluctuation  level  of  Sum,  AS,  due  both  to  statist  cal 
variability  of  the  returned  signal  (so-called  ’speckle  noise")  and  system 
noise,  be  less  than  the  variation  in  Son,  due  to  some  prescribed  variation 
(desired  precision)  in  vz,  6vz;  i.e., 

(AS)2  <  2Snin5S  •  2Suiii(9S/9vz)max^vz  *  4<qLSvz/v  ,  (32) 

where  the  subscript  q  refers  to  the  center  frequency  of  the  transmitted 
(narrowband)  pulse.  Equation  (32)  has  been  evaluated  at  Tmax  311(1  snm  has  been 
taken  as  a  i.  If  N_1  is  the  system  signal-to-noise  ratio,  then  (AS)2  is  given 


<(AS)2>  a  2(N2  +  (1/BTp) ]  ,  133 

the  factor  2  arising  because  S^  is  a  second-order  statistic,  whence  equation 
(32)  becomes 

[(1/BTp)  +  N2 J  §  [2KqL(6vz/v*)J  .  (3‘ 
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Equation  (34)  can  be  further  simplified.  Let  1  be  the  spatial  scale  of  the 
velocity  field  necessary  to  be  resolved.  A  pulse  duration  time  Tp  must  then 
be  chosen  that  is  no  longer  than 

Tp  -  2  1/c  .  (35) 

For  a  narrow-bandwidth  pulse,  one  can  write 

B  s  uq/Q  ,  (36) 

where  o>q  is  the  center  frequency  of  the  transmitted  pulse;  for  a  narrow- 
beamwidth  planar  transducer,  one  can  use  the  approximation  (Kinsler  et  a  1 . , 
1982) 


L/(2ir)  -  1/iCqU  ,  (37) 

where  u  is  the  half-angle  beamwidth,  Kq  -  wq/c.  Substituting  equations  (35) 
through  (37)  into  (32)  yields 

[Q/(2fcql )  +  N2}  3  [4ir  5vz/(uv*)]  .  (38) 

Equation  (38)  can  be  considered  the  fundamental  equation  for  setting 
sensitivity/resolution  limits  for  the  mean  Doppler  estimate  vz. 

From  figure  3,  the  model  wavenumber  spectrum  of  velocity  microstructure 
goes  as  Ek-5^.  The  velocity  sensitivity  Avz  corresponding  to  the  spatial 
scale  is  given  by 

6vz  s  E' 11/3  ,  (39) 

where,  from  figure  3,  E'  has  the  typical  range  of  values 
3  x  10~3m2,/3/sec  >  E'  >  10_3m2>/3/sec  , 
corresponding  to  strong  and  weak  turbulence  limits  ( i .  e . ,  curves  I  and  III, 
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respectively,  in  figure  3).  Substitution  of  equation  (39)  into  (38)  for  the 
case  of  zero  system  noise  (N  -  0)  yields  the  minimum  resolution  criterion 

1  £  [Quv V (8tncqC ' ) ]0 •  75  ,  (40) 

while  for  the  case  dominated  by  system  noise,  it  yields 

1  £  [N2uv*/(4irE')]3  .  (41) 

As  an  example,  for  the  sonar  parameters, 

u  ■  1°  -  1.7  x  10-2  radians, 

Q  -  15,  and 

fg  -  100  kHz  and  1  MHz, 

table  la  shows  the  minimum  spatial  resolution  for  a  signal-to-noise  ratio  of 
20  log  (N)  »  -10  dB.  A  value  of  v*  -  0.3  meter/second  and  the  weak  turbulence 
limit  of  figure  3  (E '  -  lO-3!!!2^3 1 1/f3 )  have  been  used.  Table  la  indicates  that 
over  this  frequency  range  coherent  Doppler  does  have  sufficient  spatial 
resolution  to  resolve  the  microstructure  spectrum  to  scales  on  the  order  of 
centimeters.  However,  to  achieve  this  .  ^1  resolution  for  an  isotropic 

microstructure  field,  the  transverse  dimension  of  the  scattering  volume 
(perpendicular  to  the  mean  ray  axis)  must  be  on  the  order  of  1  and,  hence,  the 
range  is  restricted  to 

R  <  l/(2u)  s  30  1  , 

with  the  range  limited  to  less  than  1  meter  for  1  on  the  order  of 
centimeters.  However,  the  range  extent  of  the  nearfield  of  a  planar 
transducer  can  be  written  (Kinsler  et  al . ,  1982)  as 

^nearfield  *  w/(8tcqu2)  s  3.25  meters  and  0.325  meters 
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Table  la.  Sensitivity  Limits:  Spatial  Resolution  of  Coherent  Doppler 


N2  -  0 

N2  -  0.1 

1  MHz 

0.44  cm 

0.60  cm 

100  kHz 

2 .5  cm 

3.1  cm 

for  the  100-kHz  and  1-MHz  cases,  respectively.  Thus,  although  the  sonar  has 
the  spatial  resolution  and  velocity  sensitivity  necessary  to  resolve 
microstructure  using  a  coherent  Doppler  estimate,  the  range  over  which  this 
occurs  is  very  restrictive.  The  upper  limit  of  l  MHz  for  the  calculation  used 
here  was  chosen  somewhat  arbitrarily  based  on  transducer  size,  power  output, 
and  range  constraints.  However,  both  equations  (40)  and  (41)  show  a  less  than 
linear  increase  in  spatial  resolution  with  increasing  frequency. 

Note  that  equations  (23)  and  (25)  establish  relationships  between  the 
statistics  of  the  received  signal  and  those  of  the  velocity  field,  namely, 
<(v'z)2>d  and  2[l  -  Rv(z)]<(v'z)2>d*  the  variance  and  the  structure  function, 
averaged  over  the  scattering  volume.  Since  this  average  is  taken  over  a 
spatial  scale  on  the  order  of  the  largest  scale  of  the  microstructure  field, 
these  relationships  allow  a  relaxation  of  the  spatial  resolution  criteria  used 
in  Doppler  estimates  of  the  instantaneous  velocity  fluctuation  from  the 
smallest  scale  of  the  variability  field  to  the  largest  scale  of  the 
variability,  provided  that  the  sensitivity  of  the  estimate  of  the  velocity 
field  statistics  is  sufficient  to  resolve  the  microstructure.  This,  in  turn, 
would  increase  the  range  capability  to  the  order  of  30  meters  for  1  -  l  meter 
and  to  300  meters  for  1  -  10  meters.  The  estimates  of  <(v'z)2>d  and 
2[l  -  Rv(z)]<(v'z)2>d  obtained  through  the  use  of  equations  (23)  and  (25)  are 
called  stochastic  Doppler  estimates. 
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From  this  viewpoint,  the  sonar  parameters  are  set  such  that  the 
scattering  volume  is  on  the  order  of  (or  larger  than)  the  characteristic 
volume  of  the  velocity  microstructure,  which  is  on  the  order  of  meters.  Note 
that,  for  1  >  1  meter,  the  mean  velocity  component  vz  corresponds  to 
oceanographic  processes  of  larger  scale  than  microstructure  and,  from  equation 
(38),  can  be  estimated  to  be  on  the  order  of  0.1  mm/sec  for  N2  -  0.1  and  for 
the  same  sonar  parameters  used  in  the  calculations  of  table  la.  Thus,  the 
variance  of  a  vz  estimate  corresponding  to  spatial  scales  greater  than  1  meter 
would  not  be  a  limiting  factor  in  the  estimation  of  parameters  associated  with 
the  microstructure  field  by  the  use  of  equations  (23)  and  (25).  A 
characteristic  parameter  that  has  arisen  in  the  stochastic  Doppler 
relationships  of  (23)  and  (25)  is  the  decorrelation  time  (namely,  tq) 
associated  with  the  variability  of  scales  on  the  order  of  the  scattering 
volume  and  smaller  given  by  equation  (24a).  The  maximum  number  of  correlation 
times  obtainable  is  limited  by  the  effect  of  transverse  flow  past  the  sonar 
and,  following  the  discussion  above,  is  easily  seen  as  given  by 

^max  -  Tmax/'t0  -  [2<(v’z)2>(|Jl''2icqL/v* 

*  (2ir)[2<  (v’z)2>(jl1/2/(uv*)  s  1.74  (42) 

for  the  sonar  parameters  used  above,  with  v*  -  0.3  meter/second  and  <(v  z)2>d 
calculated  from  equation  (39)  and  the  weak  turbulence  limit.  To  obtain  the 
precision  possible  using  the  stochastic  Doppler  approach,  the  variance  of  Qq 
will  be  used  as  the  fundamental  relationship  between  the  acoustics  and  the 
fluid  velocity  field,  with  the  along-beam  mean  component  removed  by  equation 
(18).  Accordingly,  from  equations  (23)  and  (25), 

<(^><1  “  /J>  d£d£'g(l- )g(£ ’ )  {exp[-2T2RV(z-z ' ) 3 — 1 }  , 
where  t  s  xmax  ,  with  xmax  given  by  equation  (42). 

From  equation  (29)  with  1  -  AZ  for  x  -  xmax  , 

<Qq2>d  3  V^rnax1)  -  { lvv*v>[2<v' z )2 >d ] 1 /2 / ( 2ir  1 ) >  .  (43) 
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Thus,  equation  (43)  can  be  used  as  the  smallest  <Qq2>d  to  be  resolved  in  order 
to  have  sufficient  sensitivity  to  resolve  the  smallest  scale  of  the 
microstructure  lv.  Let  (AS)2  be  the  fluctuation  level  of  (Qq)2  due  to  speckle 
noise  and  system  noise;  then,  since  equation  (43)  is  a  fourth-order  statistic, 

(AS)2  -  4(l/BTp  +  N2)  .  (44) 

The  sensitivity  can  be  further  improved  if  there  exists  some  number  of 
independent  samples  of  the  returned  signal  no  over  which  v’z  is  stationary,  in 
which  case  (AS)2  would  be  reduced  by  l/no.  Thus,  the  criterion  for  resolving 
the  smallest  scale  of  microstructure  lv  is 

(l/n0) (AS)2  <  Uv/^maxU]2  •  (45) 

One  can  write 

°0  ■  Ta/,Tmax  *  Vx^/L  -  al/L  ,  (46) 

where  xa  -  al/v*  is  the  time  period  over  which  v'z  is  statistically 
stationary.  The  parameter  a  is  associated  with  the  degree  of  anisotropy  of 
the  larger  scale  variability  generating  the  microstructure.  If  this  field  has 
physics  similar  to  that  of  internal  waves,  one  would  expect  that  a  is  a  ratio 
of  horizontal-to-vertical  scale  of  the  larger  scale  field  and  that  a  s  N/f  • 
100,  where  N  is  the  buoyancy  frequency  and  f  is  the  coriolis  frequency.  (Note 
that  a  can  be  considered  as  the  ratio  of  the  horizontal  extent  to  the  vertical 
extent  of  a  patch  of  microstructure.)  Substituting  equations  (30),  (35)-(37), 
(44),  and  (46)  into  equation  (45)  yields 

lv  ^  {[N2  +  Q/(2iCql)][64Tr3l<(v'z)2>(j]/[K:qau3(v*)2]}1/2  .  (47) 

For  the  same  sonar  parameters  as  used  for  the  coherent  Doppler  case  (i.e., 

Q  »  15,  v  -  1°,  and  v*  -  0.3)  but  with  1  -  1  meter,  the  minimum  resolvable  lv 
for  the  l-MHz  and  100-kHz  cases  has  been  calculated  and  is  shown  in  table  lb. 
These  results  suggest  that  there  is  sufficient  resolution  to  obtain  lv, 
although  it  is  somewhat  marginal  at  100  kHz  for  the  N2  -  0.1  case.  Note  that 
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equation  (47)  depends  only  on  the  square  root  of  a,  the  degree  of  anisotropy 
of  the  microstructure  patch.  A  factor  of  10  reduction  in  a  yields  a  factor  of 
3  reduction  in  sensitivity  with  acceptable  values  in  all  cases  except  the 
100-kHz,  N2  -  0.1  case.  Also,  note  that  equation  (47)  shows  no  dependence  on 
1  for  the  noise-free  case  —  a  rather  surprising  result.  Significant 
improvement  in  sensitivity  can  occur  by  decreasing  the  mean  transverse  current 
(linear  dependence).  The  value  of  v*  -  0.3  meter/second  is  an  upper  limit 
estimate.  Finally,  although  equation  (47)  shows  an  increase  in  sensitivity 
with  increasing  beamwidth,  the  approximation  used  in  equations  (40)  and  (42) 
involved  Tmax  >  1,  which  from  equation  (42)  shows  a  decrease  in  xmax  with 
increasing  beamwidth. 

Table  lb.  Sensitivity  Limits:  Spatial  Resolution  of  Stochastic  Doppler 


N2  -  0 

N2  -  0.1 

1  MHz 

4.26  cm 

1.1  cm 

100  kHz 

13.0  cm 

3.5  cm 
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V.  SUMMARY  AND  CONCLUSIONS 


This  report  has  examined  the  returned  signal  statistics  associated  with 
acoustic  scattering  from  a  large  collection  of  scatterers  embedded  in  a  model 
fluid  velocity  field  appropriate  for  oceanographic  application.  The  fluid 
velocity  field  is  three-dimensional  and  contains  a  microstructure  or 
"turbulent"  component.  The  underlying  acoustic  model  is  quite  simple  —  a 
monostatic,  planar  transducer  projects  very  high  frequency  sound  (>100  kHz) 
into  a  very  large  number  of  identical  scatterers  that  follow  the  fluid 
perfectly.  Identical  pulses  are  transmitted  at  some  constant  repetition  rate 
Tr  and  then  gated  at  some  delay  time  corresponding  to  some  prescribed  range  of 
interest.  The  propagation  is  assumed  to  be  in  the  geometric  acoustics  regime 
and  only  the  farfield  case  is  considered.  Ambiguous  scattering  volumes 
associated  with  a  repetitive  set  of  transmitted  pulse  trains  are  ignored. 

(This  problem  is  addressable  either  by  restricting  the  repetition  rate  Tr  to 
be  much  greater  than  the  travel  time  to  a  range  where  the  reverberation  level 
is  below  some  prescribed  amount  or  by  using  coded  pulses.) 

The  components  of  velocity  that  affect  the  Doppler  statistics  of  the 
returned  signal  are: 

1.  The  velocity  component  in  the  direction  of  the  beam  pattern  axis  and 
on  the  scale  of  the  scattering  volume  and  larger  vz; 

2.  The  velocity  component  transverse  to  the  direction  of  the  beam  pattern 

axis ; 


3.  The  velocity  shear  field  in  the  direction  of  the  beam  pattern  axis; 

4.  The  velocity  components  in  the  direction  of  the  beam  pattern  axis  and 
on  the  scale  of  the  scattering  volume  and  smaller. 

Since  these  effects  have  been  addressed  by  many  authors  but  in  a  variety  of 
contexts,  a  detailed  derivation  of  the  relationship  of  the  returned  signal 
statistics  in  terms  of  these  velocity  components  is  presented  in  the 


appendix;  the  derivation  can  be  summarized  through  equations  (13),  (15),  and 
(16). 


The  along-beam  velocity  on  the  scale  of  the  scattering  volume  and  larger 
produces  a  Doppler  shift  or,  equivalently,  a  time  delay  in  the  pulse-to-pulse 
correlation  function  S^.  The  transverse  component  of  flow  produces  Doppler 
broadening  or,  equivalently,  a  decorrelation  in  Skq.  The  decorrelation  time 
Tmax  associated  with  the  transverse  component  of  flow  is  shown  to  be  on  the 
order  of  L/v*,  where  L  is  a  characteristic  dimension  of  the  sonar  in  the 
direction  of  the  mean  transverse  flow  velocity  v*.  This  decorrelation  time 
sets  a  maximum  time  over  which  a  moored  Doppler  sonar  can  receive  coherent 
information  for  the  model  scattering  field  presented  in  this  work. 

Decorrelation  effects  or,  equivalently,  Doppler  spread  effects,  caused  by 
velocity  variability  (microstructure)  on  scales  smaller  than  the  scattering 
volume  show  up  via  the  statistics  of  Qi ,  equation  (16),  and  have  a 
characteristic  decorrelation  time  tq.  given  by  equation  (24b).  Calculations 
in  section  IV,  i.e.,  equation  (42),  using  typical  oceanographic  values  and 
sonar  parameters  show  that  xmax/xo  s  1.7,  and  that  this  ratio  is  independent 
of  frequency  but  dependent  on  the  beamwidth  and  the  ratio  of  the  micro- 
structure  velocity  to  the  transverse  velocity  v*.  These  results  indicate  that 
there  is  enough  time  of  coherence  to  obtain  the  variance  of  the  velocity 
microstructure. 

In  section  IV  a  comparison  was  made  between  a  coherent  Doppler  estimate 
and  an  estimate  based  on  the  statistical  relationships,  equations  (23)  and 
(25),  called  stochastic  Doppler.  The  coherent  Doppler  estimate  from  frequency 
ranges  of  100  kHz  to  1  MHz  has  sufficient  sensitivity  to  resolve  the  smallest 
scales  of  ocean  velocity  microstructure  but  is  very  restrictive  in  range.  On 
the  other  hand,  over  the  same  frequency  range,  the  stochastic  Doppler 
technique  also  appears  to  have  the  sensitivity  to  resolve  the  smallest  scales 
of  ocean  velocity  microstructure  with  much  less  restriction  on  range. 


30 


VI.  BIBLIOGRAPHY 


Albright,  R.J.  (1976),  "Relationship  of  Doppler  Ultrasonic  Scattered  Signal 
Characteristics  to  Flow  and  Beam  Patterns,"  Journal  of  the  Acoustical 
Society  of  America,  vol  59,  no. 4. 

Brown,  E.H.  (1974),  "Turbulent  Spectral  Broadening  of  Backscattered  Acoustic 
Pulses,"  Journal  of  the  Acoustical  Society  of  America,  vol  56,  no.  5. 

Brown,  E.H. ,  and  S.C.  Clifford  (1973),  "Spectral  Broadening  of  an  Acoustic 
Pulse  Propagating  Through  Turbulence,"  Journal  of  the  Acoustical  Society 
of  America,  vol  54,  no.  1. 

Edwards,  R.V.,  J.C.  Angus,  and  J.W.  Dunning,  Jr.  (1971),  "Spectral  Analysis 

of  the  Signal  from  the  Laser  Doppler  Flowmeter:  Time  Independent  System," 
Journal  of  Applied  Physics,  vol  24,  pp  837-850. 

Farmer,  D.M.,  and  G.D.R.  Crawford  (1983),  "Measurements  of  Acoustic 

Correlation  in  the  Ocean  with  a  High  Frequency  Echosounder ,"  Nature ,  vol 
36,  no.  580,  pp  1-3. 

Flatte,  S.,  R.  Dashen,  V.  Munk,  K.  Watson,  and  F.  Zachariasen  (1979),  Sound 

Transmission  Through  a  Fluctuating  Ocean,  Cambridge  University  Press,  New 
York. 

Gargett,  A.E.,  P.J.  Hendricks,  T.B.  Sanford,  T.R.  Osborn,  and  A. J. Williams 
(1981),  "A  Composite  Spectrum  of  Vertical  Shear  in  the  Upper  Ocean," 
Journal  of  Physical  Oceanography,  vol  11,  pp  1258-1271. 

Green,  S.  (1964),  "Spectral  Broadening  of  Acoustic  Reverberation  in 

Doppler-Shift  Fluid  Flowmeters,”  Journal  of  the  Acoustical  Society  of 
America,  vol  36,  no.  7. 

Hinze,  J.O.  (1959),  Turbulence,  McGraw-Hill,  New  York. 

Joyce,  T.,  D.  Bitterman,  Jr.,  and  K.E.  Prada  (1982),  "Shipboard  Acoustic 

Profiling  of  Upper  Ocean  Currents,"  Deep  Sea  Research,  vol  29,  no.  7A,  pp 
903-913 

Kinsler,  L.E.,  A.R.  Frey,  A.B.  Coppens,  and  J.V.  Sanders  (1982),  Fundamentals 
of  Acoustics,  John  Wiley  and  Sons,  New  York. 

Lamb,  H.,  1932,  Hydrodynamics,  Cambridge  University  Press,  London. 

L'Hermitte,  R.  (1983),  "Pulse  Coherent  Doppler  Sonar  for  the  Measurement  of 
Water  Current  and  Turbulence,"  Proceedings  of  the  Acoustic  Profiling 
Symposium,  November  1983,  IEEE  Current  Measurement  Tech  Committee,  NOAA. 

Mathews,  T.  and  R.B.  Hicks  (1981),  "Proceedings  of  the  International  Symposium 
on  Acoustic  Remote  Sensing  of  the  Atmosphere  and  Oceans,"  University  of 
Calgary,  Alberta,  Canada. 


31 


Newhouse,  V.L.,  E.C.  Furgason,  G.F.  Johnson,  and  D.A.  Wolf  (1980),  "The 
Dependence  of  Ultrasound  Doppler  Bandwidth  on  Beam  Geometry,"  IEEE 
Transactions  on  Sonics  and  Ultrasonics,  vol  27,  no.  2. 

Papoulis,  A.  (1965),  Probability,  Random  Variables,  and  Stochastic  Processes, 
McGraw-Hill,  New  York. 

Pinkel,  R.  (1979),  "Acoustic  Doppler  Techniques,"  MPL-U-86/76,  Marine  Physical 
Laboratory,  San  Diego,  CA. 

Pinkel,  R.  (1980),  "On  the  Use  of  Doppler  Sonar  for  Internal  Wave 

Measurements,"  MPL-U-63/79,  Marine  Physical  Laboratory,  San  Diego,  CA. 

Pinkel,  R.  (1981),  ”0n  the  Use  of  Doppler  Sonar  for  Internal  Wave 
Measurements,"  Deep  Sea  Research,  vol  28A,  no.  3,  pp  269-280. 

Rasmussen,  R.A.,  and  N.E.  Head  (1978),  "Remote  Detection  of  Turbulence  from 

Observations  of  Reverberation  Spectra,"  Journal  of  the  Acoustical  Society 
of  America,  vol  63,  no.  1. 

Regier,  L.A.  (1982),  "Mesoscale  Current  Fields  Observed  with  a  Shipboard 

Acoustic  Profiling  Current  Meter,"  Journal  of  Physical  Oceanography,  vol 
12,  no.  8,  pp  880-886. 

Rummler,  W.D.  (1968),  "Introduction  of  a  New  Estimator  for  Velocity  Spectral 
Parameters,"  Technical  Meow  Mm-68-4145-5,  Bell  Telephone  Laboratories, 
Whippany,  NJ. 

Sullivan,  E.J.,  and  K.A.  Kemp  (1979),  "The  Significance  of  Beam  Pattern  Side 
Lobes  in  Fluid  Flow  Measurement,”  Journal  of  the  Acoustical  Society  of 
America,  vol  66,  no.  5. 

Tatarski,  V.I.  (1971),  "The  Effects  of  Turbulent  Atmosphere  on  Wave 

Propagation,"  Israel  Program  for  Scientific  Translation,  Jerusalem. 


32 


APPENDIX 

DERIVATION  OF  PULSE-TO-PULSE  CORRELATION  FUNCTION 

As  discussed  in  the  body  of  this  report,  the  first-order  statistics  of 
the  received  signal  ^(t')  are  zero,  as  would  be  expected.  For  the  second- 
order  statistics, 

®nm  ”  ^n^m  • 

Then,  the  pulse-to-pulse  gated  autocorrelation  function,  equation  (11),  is 
given  by 

Snm  “  [h2( t  )®nmlTi  <  ( A— 1 ) 

wi  th 

®nm  -  a(Rn)a(Rm)  S  dr’  /  dr'^n’ )b2(n" )  E  W*1 
Vn  Vm  k,l 

exp{i  (u^i ) t '  +2 i [«k(Rn-r ’  )-«l (Rm-r" ) ]/c} 

p(r',t'  +  (n  -  l)Tr  +  Tjj]  p[r" ,  t  ’  +  (m  -  l)Tr  +  xm)  , 

and  the  scattering  volumes  Vn  and  Vm  are  limited  in  radial  extent  by 

rj  £  r'  -  Rn  S  r2  for  Vn  , 

and 

U  S  r'  -  Rn  S  r2  for  Vm  , 

with 


rl  -  c( t ’ 

-  Tp/2)/2  , 

(A-2a) 

f2  -  c  ( t  ’ 

-  Tp/2 )/2  . 

(A-2b) 

A-l 

Since  the  initial  position  of  each  scatterer  Xj  is  taken  as  equiprobable,  the 
density  function  p,  being  a  sum  of  delta  functions,  can  be  considered  from  the 
definition  of  the  Lagrangian  position  vector  r,  i.e.,  from 


r(Xj ,  t)  «  Xj  +•  dj  , 

as  a  Poisson  process  with  respect  to  Xj  (Albright,  1976).  Note,  also,  that 
the  scatterer  displacement  is  assumed  to  be  small  enough  that  the  explicit 
time  dependence  in  p(r,t)  can  be  neglected  over  the  pulse  duration  time  Tp  and 
the  range  gate  time  Tf.  Since  the  scatterer  initial  position  Xj  can  be 
described  by  a  Poisson  process,  it  is  straightforward  to  show  that  to  order 
1/N,  where  N  is  the  average  number  of  scatterers  in  the  scattering  volume, 

Snm  “  [h2 ( t  ’  )©njnJT|  a  <U2(t '  )®nm]Tj>X  -  [h2(t '  )<9nm>x]Tr  .  (A-3 ) 


where 


<0nm>X  -  a(Rn)a(Rm)  X  dr’  X  dr’b2^'  )b2(n" )  E  W*1 

Vn  Vm  k.l 

exp  { i  i )  t  ’ +2  i  (wfe  ( Rn- r ' ) -« i  ( Ra,- r  ” )  ] /c } 

<p[r',t'  +  (n  -  l)Tr  +  xn)  p(r",t'  +  (m  -  l)Tr  +  Tm]>x 

-  a(Rn)a(Rm)  X  dr'  X  dr"b2(n  )b2(n”)  E  W'\ 

Vn  Vm  k.l 

exp{i(«jc-mi)t,+2i[uk(Rn-r’ )-«l(Rm-r")]/c}  po$(*'  -  r"  +  s). 

(A-4) 

Note  that  the  second-order  statistic  for  a  Poisson  process  (Papoulis,  1965) 
has  been  used;  namely, 


A-2 


< P [ r ' , t ’  +  (n  -  1 )Tr  +  xnJ  ptr",!’  +  (m  -  l)Tr  +  xm]>x 


s  <p[r',(n  -  l)Tr  +  xn]  p[r",(m  -  l)Tr  +  xm]>x 
-  PO2  +  POs(r'  “  r"  +  s)  -  ( A— 5 ) 

where 


s  «  d[X, (m  -  l)Tr  +  Tqj]  -  d[X, (n  -  1)T^  +  x^] 

is  the  displacement  of  a  fluid  element,  initially  at  X,  during  the  time 
interval  between  the  scattering  of  the  nth  and  mth  pulse.  Again,  since  it  has 
been  assumed  that  the  pulse  duration  Tp  and  the  gate  function  Tj  are  kept 
sufficiently  small,  the  explicit  time  dependence  of  p  on  t’  can  be  neglected 
in  equation  (A-4).  However,  since  Tr  »  Tp  +  Tj  ,  the  dependence  on 
pulse-to-pulse  time  delay  cannot  be  ignored.  Note  that  the  first  term  on  the 
right-hand  side  of  equation  'A-5)  does  not  contribute  to  (A-4)  by  virtue  of 
equation  (10).  (The  product  of  two  zero-mean  value  terms  is  zero.) 

Consider  the  effect  of  the  delta  function  on  the  beam  pattern  functions 
of  equation  (A-4).  For  nonzero  contributions,  the  delta  function  in  (A-4) 
requires  that 

r"  -  r’  +  s  . 


Then, 

n"  -  (r1  +  s) / 1 r *  +  sj  a  n'  ,  (A-6) 

since  the  displacements  of  the  scatterer  over  the  time  intervals  between 
pulses  are  such  that 

|s|  «  |f|  . 

For  the  effect  of  the  delta  function  on  the  radial  spatial  variables  of 

A-3 


equation  (A-4),  note  that  contributions  can  occur  only  if  the  double  integral 
overlaps  and,  hence,  that 

| r*  -  r"  +  s|  s  cTp/2  . 

The  condition  for  overlap  of  the  radial  component  of  the  spatial  integrals  of 
equation  (A-4)  can  be  written 

Rm  "  Rn  +  cAx/2  -  R  +  cAx/2  ,  ( A— 7 ) 

where  R  is  now  identified  as  the  range  to  the  nth  scattering  volume,  used  as 
a  reference.  The  time  delay  Ax  will  be  used  to  estimate  a  mean  Doppler 
shift.  Integrating  over  the  delta  function  in  equation  (A-4)  using  (A-6) 
results  in 

<®nm>X  -  R2*2(R)po  X  drb4(n)  l  1 

V0  k,l 

exp{i(u|c-ui)(t'+2(R-r)/c]-i(i»i(Ax-2Sr/c)}  ,  ( A-8 ) 

where  Vq,  the  scattering  volume,  constrains  the  range  interval  by 
rx  £  r  -  R  3  r2  , 

with  r1(  r2  defined  by  equations  (A-2a)  and  (A-2b),  respectively.  The  radial 
displacement  Sr  is  given  by 

Sr  -  n  •  s.  ( A— 9 ) 

Note  that  in  the  dummy  spatial  variable  r  of  integration  in  equation  (A-8)  the 
prime  notation  has  been  dropped. 

The  beam  pattern  function  can  be  written  in  the  farfield  as  the  Fourier 
transform  of  the  aperture  function  w(x);  namely, 

b(n)  -  S  dxw(x)  exp( iwqn-x/c)  ,  (A-10) 
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where  attention  is  confined  to  a  planar  sonar  and  x  is  taken  in  the  x,y  plane, 
i.e.,  x  =  (x,y,0).  (Note  that  the  effect  of  finite  frequency  bandwidth 
dependence  of  the  beam  pattern  has  been  neglected;  this  effect  can  easily  be 
taken  into  account  for  the  narrowband  pulses  in  the  analysis  by  replacing  b(n) 
with  bfc(n),  where  k  indicates  the  pulse  Fourier  series  component  number.  The 
frequency  uq  is  the  center  frequency  of  transmission.) 

Note  that  for  a  planar  sonar 
n  •  x  =  nh  •  x  , 

where  %  is  the  vector  component  of  the  radial  unit  vector  n  in  the  x,y  plane, 
i.e., 

nh  -  (nx,  ny ,  0) . 


Then, 

b(n)  -  b(nh) .  (A-ll) 

Substitution  of  equation  (A-ll)  into  (A-8)  yields 

<9nm>X  ”  R2<*2(R)P0  J*  drb4(nh)  Z  W*1 

V0  k,l 

exp{i(ujc-ui)[t'+2(R-r]/c-iui(Ax-2Sr/c)}  .  ( A-12 ) 

For  a  narrowbeam  sonar, 
nz  -  cos(u)  -  1  , 

since  the  polar  angle  u,  defined  with  respect  to  the  mean  ray  axis  (taken  as 
the  z-direction) ,  is  small,  i.e.,  u  «  1.  Then,  the  radial  spatial  variable  of 
integration  in  equation  (A-12)  approximates  to 
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r  s  r  cos(u)  -  z  , 


(A— 13 ) 


over  the  limits  given  by  (A— 9) : 


ri  ^  z  -  R  ^  f2 


( A— 14 ) 


where  r^  r£  are  given  by  equations  (A-2a)  and  (A-2b),  respectively.  By 
definition,  the  spatial  variables  of  integration  of  equation  (A-12)  in  the  x,y 
plane  can  be  written  as 

x  •  rnjj  s  znjj  (A-15 ) 

by  virtue  of  (A-13).  Using  equations  (A-13)  and  (A-15),  the  spatial  volume 
element  of  (A— 9 )  can  be  written  as 

d3x  s  z2dzdnh  3  R2dzdi%  , 

which  results  in  (A-12)  becoming 

<®nm>X  •  R2a2(R)po  /  di%b4(i%) 

R+c(  t ' +Tp/2)/2 

E  W*1  X  dz  exp(i(»ic-«i)[t  +2(R-z)/c]-iMi(Ax-2Sr/c)} .  (A-16) 
1 ,k  R+c(t’-Tp/2)/2 


Let 


f  1  (z )  -  exp(2i«iSr/c)  , 


where 


-cTp/4  <  z  -  R  -  ct'/2  <  cTp/4  .  (A— 17 ) 

To  simplify  equation  (A-16),  fi(z)  is  expanded  over  (A-17)  in  the  Fourier 
series 
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( A— 18 ) 


f  1  ( z )  -  exp(2i«iSr/c)  -  Z  aij  exp[ i (z-R-ct ' /2)tc j ]  , 

j 

where 

tcj  -  2ir j/ (cTp/2)  , 
and 

R+ct ' /2+cTp/4 

ajj  *  l/(cTp/2)  J  dz  exp(2iuiSr/c)  exp[-  i  (z-R-ct '/2)icj]  ,  ( A— 19 ) 

R+ct ' /2-cTp/4 

which,  upon  substitution  into  equation  (A-16)  and  after  some  straightforward 
algebra,  yields 

<0nm>X  *  R2a2(R)p0  J  dnhb4(i%)(cTp/2)  Z  vi+jY*i 

I.  j 

exp(-iuidT)ai t j  ,  (A-20) 

and  then  substituting  equation  (A-20)  into  (A— 3 )  yields 

Snm  “  R2a2(R)P0  f  dnhb4(nh)(cTp/2)  Z  tl+jt*l 

l.j 

exp(-iuiA*c)[h2(t’  )ai,  j]Tp  .  (A-21) 

Equation  (A-20)  shows  that,  in  general,  individual  Fourier  components  of  the 
transmitted  pulse  mix,  i.e.,  tl+jt*!-  contribute  for  m  #  0  to  the  correlation 
function  S^.  Note  that  for  homogeneous  Gaussian  displacement  statistics, 

<Sr (z ’ )Sr(z" )>d  -  F(z”  -  z)  , 

<al ,  j>d  ■  0  for  j  *  0  ; 
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thus,  for  <Snffl>,  no  mixing  of  pulse  Fourier  components  occurs.  So,  if  one 
chooses  a  scattering  volume  sufficiently  large  to  include  all  the  scales  of 
variability,  one  expects  that 

(al, j>Tp  -  < a 1 , j >d  -  0  for  j  n  0  .  (A- 

If  one  assumes  condition  (A— 22 )  in  equation  (A-21), 

Snm  =  R2a2(R)p0  S  dnhb4(nh)(cTp/2)ifq<j»*i 

exp(-iuiAT)[h2(t' )ai)0]Tp.  (A- 

It  should  be  noted  that  it  is  straightforward  to  include  the  j  *  0  terms  in 
equation  (A-23)  and  in  the  analysis  of  section  III.  Note  that  the  last  term 
of  (A-23)  can  be  written  for  h  -  1  as 

Tp/2  R+ct ' /2+cTp/ 4 

[  h2  ( t '  )aii0]ip  -  (a^oJTp  X  h2dt'  X  dz  exp(2ia»iSr/c) 

-Tp/2  LR+ct'/2-cTp/4 

-  X  g(C)dC  exp[2iuiSr(z)/c]  ,  (A- 

where  g(C)  is  the  triangle  function: 

g(S)  -  1  -  5  for  0  S  £  *  1  , 

-  £  -  1  for  -l  ^  i  £  0  , 

-  0  otherwise  , 


and  where  C  -  (z  -  R)/Az,  with  Az  -  cTp/2  .  Thus, 


Smn  -  R2a2(R)Tpp0  I  S  dnhb4(nh)<m*  1  S  g(5)dC 


exp{-ia>i[Ax-2Sr(z)/c]}  .  (A-25) 

The  form  of  equation  (A-25)  suggests  decomposing  the  displacement  field  into  a 
mean  and  fluctuating  component  as  follows 

s  -  s  +  s'  ,  (A-26) 

where  the  overbar  notation  is  used  to  indicate  the  weighted  spatial  average: 

F(z]  -  J*  g(£)f(z) . 

Thus,  by  definition, 

s  -  d[X, (m  -  1 ) Tr  +  xm]  -  3[X, (n  -  l)Tr  +  xn]  , 

■  v(m  -  n)Tr  ,  (A-27 ) 

where  equation  (A-7)  has  been  used  and,  hence,  Tm  -  xn  +  Ax,  Ax  «  xn,  xm  <  Tr, 
and  where  the  time  interval  (m  -  n)Tr  has  been  selected  to  be  small  enough 
that  the  mean  velocity  v  does  not  change  significantly  over  (m  -  n)Tr.  The 
fluctuating  displacement  s'  can  be  written  as 

s'  -  d  [X, (m  -  l)Ty  +  Xn]  -  d'[X,(n  -  l)Tr  +  x^]  ( A— 28 ) 

Note  that  for  turbulence- like  fluctuations,  if  (m  -  n)Tr  <  x\  ,  where  x\  is 
the  temporal  microscale,  the  frozen  field  hypothesis  (Hinze,  1959)  can  be 
invoked  and  s'  can  also  be  written  in  the  form 

s'  -  (m  -  n)Trv'(X)  ,  (A-29) 

where  v'  is  the  fluctuating  velocity  defined  analogously  to  equation  (A-27). 
Substituting  equations  (A-27)  and  (A-29),  using  (A-10),  into  equation  (A-25) 
yields 
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Snjn  -  R2a2(R)Tppo  E  1  X  di%b4(nh)  exp[ i2o>ivh(m-n)Tr -i%/c] 

1 

exp{iuj[2v2(m-n)Tr/c-Ax]}  X  d£g(£)  exp{iu[2v'z(m-n)Tr/c]} ,  (A-30) 

where  the  overbar  and  prime  notation  (for  mean  and  fluctuating)  have  been 
combined  with  the  vector  component  notation  subscript  h  (indicating  along  the 
x/y  plane)  and  subscript  z  (indicating  the  direction  of  the  mean  ray  axis). 

Vote  that  sh  -  vh(m-n)Tr  »  s'h-  (It  is  straightforward  to  include  the 
term  s'h  ~  v'h(m  -  n)Tr  with  some  additional  algebra.  However,  most  realistic 
oceanographic  cases  will  allow  this  generalization  since  the  oceanographic 
velocity  variability  is  typically  "red"  in  wavenumber  spectra.)  To  normalize 
equation  (A-30),  take  S^/Sna  ,  which  amounts  to  setting 

x  .  -«r  ,  • 

R2a2(R)Tppo  -  1  . 
and 

E  V1H»*1  -  1  . 

1 

After  redefining  Sqq  as  the  normalized  pulse-to-pulse  gated  correlation 
function  wi+h  Snn  »  1  from  equation  (A-30),  the  above  normalization  then  yields 

Snm  -  £  tit* 1  X  dnhb4(nh)  exp[i2uivh(m-n)Tr*nh/c] 

1 

exp{i«i(2vz(m-n)Tr/c-bx]}  X  dCg(C)  exp{i»[2v"z(m-n)Tr/c]}  .  (A-31) 
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